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Abstract 



This is a review of the authors' recent results on an integrable structure 
of the melting crystal model with external potentials. The partition func- 
tion of this model is a sum over all plane partitions (3D Young diagrams). 
By the method of transfer matrices, this sum turns into a sum over ordi- 
^ , nary partitions (Young diagrams), which may be thought of as a model 

f^ ' of q -deformed random partitions. This model can be further translated 

r~^ I to the language of a complex fermion system. A fermionic realization of 

^\ . the quantum torus Lie algebra is shown to underlie therein. With the 

^^ ' aid of hidden symmetry of this Lie algebra, the partition function of the 

melting crystal model turns out to coincide, up to a simple factor, with a 
(•—^ ■ tau function of the ID Toda hierarchy. Some related issues on 4D and 5D 

QQ ' supersymmetric Yang-Mills theories, topological strings and the 2D Toda 

f^ , hierarchy are briefly discussed. 
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Figure 1: Melting crystal corner 



1 Introduction 



The melting crystal model is a model of statistical mechanics that describes a 
melting corner of a semi- infinite crystal (Figured]) ■ The crystal is made of unit 
cubes, which are initially placed at regular positions and fills the positive octant 
x,y, z > of the three dimensional Euclidean space. As the crystal melts, a 
finite number of cubes are removed from the corner. The present model excludes 
such crystals that have "overhangs" viewed from the (1,1,1) direction. In other 
words, the complement of the crystal in the positive octant is assumed to be 
a 3D analogue of Young diagrams (FigurelJI). Since 3D Young diagrams are 
represented by "plane partitions" , the melting crystal model is also referred to 
as a model of "random plane partitions" . 

Though combinatorics of plane partitions has a rather long history T] , Ok- 
ounkov and Rcshetikhin ^ proposed an entirely new approach in the course of 
their study on a kind of stochastic process of random partitions (the Schur pro- 
cess). Their approach was based on "diagonal slices" of 3D Young diagrams and 
"transfer matrices" between those slices. As a byproduct, they could re-derive 
a classical result of MacMahon [T] on the generating function of the numbers 
of plane partitions. Actually, this generating function is nothing but the par- 
tition function of the aforementioned melting crystal model. The method of 
Okounkov and Reshetikhin was soon generalized [3] to deal with the topological 
vertex [U [5] of A-model topological strings on toric Calabi-Yau threefolds. 

The melting crystal model is also closely related to supersymmetric gauge 
theories. Namely, with slightest modification, the partition function can be 
interpreted as the instanton sum of 5D JV — I supersymmetric (SUSY) U{1) 
Yang-Mills theory on partially compactified space-time R"' x S'-'^ [6]. This in- 




Figure 2: 3D Young diagram as complement of crystal corner 



Stanton sum is a 5D analogue of Nekrasov's instanton sum for 4D Af — 2 SUSY 
gauge theories [71 \8\ . The 4D instanton sum is a statistical sum over ordinary 
partitions (or "colored" partitions in the case of SU{N) theory), hence a model 
of random partitions. Nekrasov and Okounkov [9] used such models of random 
partitions to re-derive the Seiberg-Witten solutions [10] of 4D A/" = 2 SUSY 
gauge theories. Actually, by the aforementioned method of transfer matrices, 
the statistical sum over plane partitions can be reorganized to a sum over par- 
titions. This is a kind of g-deformations of 4D instanton sums. A 5D analogue 
of the Seiberg-Witten solution can be derived from this g-deformed instanton 
sum [iln]. 

In this paper, we review our recent results [12] on an integrable structure 
of the melting crystal model (and the 5D U{1) instanton sum) with external 
potentials. The partition function Zp{t) of this model is a function of the cou- 
pling constants t = (ti, ^2, . . .) of the external potentials. A main conclusion of 
these results is that Zp{t) is, up to a simple factor, a tau function of the ID 
Toda hierarchy, in other words, a tau function Tp{t, t) of the 2D Toda hierarchy 
[l3] that depends only on the difference t — i oi the two sets t, i of time vari- 
ables. To derive this conclusion, we first rewrite Zp{t) in terms of a complex 
fermion system. In the case of 4D instanton sum, such a fermionic representa- 
tion was proposed by Nekrasov et al. [TH [9] . In the present case, we can use 
the aforementioned transfer matrices [2] to construct a fermionic representation. 
This fermionic representation, however, does not take the form of a standard 
fermionic representation of the (ID or 2D) Toda hierarchy [15l[T6]. To resolve 
this problem, we derive a set of algebraic relations (referred to as "shift sym- 
metry") satisfied by the transfer matrices and a set of fermion bilinear forms. 



(Actually, these fermion bilinear forms turn out to give a realization of "quan- 
tum torus Lie algebra".) These algebraic relations enable us to rewrite the 
fermionic representation of Zp{t) to the standard form of Toda tau functions. 

In the 4D case, a similar partition function with external potentials has been 
studied by Marshakov and Nekrasov [TTlUHj- According to their results, the ID 
Toda hierarchy is also a relevant integrable structure therein. Unfortunately, 
our method developed for the 5D case relies heavily on the structure of quantum 
torus Lie algebra, which ceases to exist in the 4D setup. We shall return to this 
issue, along with some other issue, in the end of this paper. 

This paper is organized as follows. Section 2 is a brief review of the melt- 
ing crystal model and its mathematical background. Section 3 presents the 
fermionic formula of the partition function. The method of transfer matrices is 
reviewed in detail. Section 4 deals with the quantum torus Lie algebra and its 
shift symmetries. In Section 5, we use this symmetry to rewrite the fermionic 
representation of the partition function to the standard form as a Toda tau 
function. Section 6 is devoted to concluding remarks. 

2 Melting crystal model 

2.1 Young diagrams and partitions 

Let us recall [19^ that an ordinary 2D Young diagram is represented by an 
integer partition, namely, a sequence 

A = (Ai,A2,...), Ai>A2>---, 

of nonincreasing integers A^ e Z>o with only a finite number of A,;'s being 
nonzero. Ai is the length of i-th row of the Young diagram viewed as a collection 
of unit squares. We shall always identify such a partition A with a Young 
diagram. The total area of the diagram is given by the degree 



|A|=E^« 



of the partition. 

It was shown by Euler that the generating function of the number p{N) of 
partitions A of degree N has an infinite product formula: 



J2piN)q''^l[il~q'T\ (2.1) 



N=0 n=l 



where q is assumed to be in the range < g < 1. One can interpret this 
generating function as the partition function of a model of statistical mechanics. 



^2D-£p(iV)<Z^ = 5]9l- 



.|A| 
JV=0 



in which each partition A is assigned an energy proportional to |A|, and q is 
related to the temperature T as q ^ g-const./T 



2.2 3D Young diagrams and plane partitions 

A 3D Young diagram can be represented by a "plane partition" , namely, a 2D 
array 

/ TTll 7ri2 • • 
7r21 7r22 ■ • 



of nonnegative integers Tr^ £ Z>o such that 



V 



-Kij is the height of the stack of cubes placed at the (i, j)-th position of the 
plane. We shall identify such a plane partition with the corresponding 3D 
Young diagram. The total volume of the 3D Young diagram is given by 

As an analogue of p(iV), one can consider the number pp(Af) of plane par- 
titions TT with |7r| — N. The generating function of these numbers was studied 
by MacMahon pj and shown to be given, again, by an infinite product: 



N=0 



pp{N)q'' ^ 1[{1 - q")-\ (2.2) 



The right hand side is now called the MacMahon function. In statistical me- 
chanics, this generating function becomes the partition function 



^3D=£pp(iV)g'^=E'?'^ 



N=0 n 

of a canonical ensemble of plane partitions, in which each plane partition tt has 
an energy proportional to the volume |7r|. 

We shall deform this simplest model by external potentials. To this end, we 
have to introduce the notion of "diagonal slices" of a plane partition. 

2.3 Diagonal slices of 3D Young diagrams 

Given a plane partition tt = {T^ij)f°j^i, the partition 

{^^i^^+m.)^l if m > 



Tr(m) — -, , N^ .„ „ 

is called the m-th diagonal slice of tt. These partitions {TT{m)}^^_^ represent 
a sequence of 2D Young diagrams that are literally obtained by slicing the 3D 
Young diagrams (Figures]). 





(a) 



(b) 



Figure 3: Diagonal slices (b) of plane partition (a) 



The diagonal slices are not arbitrary but satisfy the condition [31 [3] 

< 7r(-2) ^ 7r(-l) ^ 7r(0) >- 7r(l) ^ 7r(2) > , (2.3) 

where 'V denotes interlacing relation, namely, 



A = (Al, A2, ...) y ^i^ (M1:M2, ■ ■ •) 



dcf 



Al > Ml > A2 > /i2 > 



Because of these interlacing relations, a pair (T, T') of semi-standard tableaux 
is obtained on the main diagonal slice A = 7r(0) by putting "m + 1" in boxes of 
the skew diagram 7r(±m)/7r(±(m + 1)). 

By this mapping tt i~^ (T, T'), the partition function Z3D of the plane par- 
titions can be converted to a triple sum over the tableau T, T' and their shape 
A: 



Z' 



3D = E E ^V, 

X T,T':shapcA 



(2.4) 



where 



= TT 9^™+^/^^''' 



(—m)/7T( — m—l) 



m=0 



.,T' 



TT (ni+l/2)|7r(m)/;r(m+l) 



m=0 



By the well known combinatorial definition of the Schur functions J19j . the 
partial sum over the semi-standard tableaux turn out to be a special value of 



the Schur functions: 

E 9^= E 'Z^'--A(g''), (2.5) 

T:shapc A T' :shapc A 

where 

Thus the partition function can be eventually rewritten as 

^3D-E^^(9')'- (2-6) 

A 

Let us note that the special value of the Schur functions has the so called Hook 
formula [19] 

SAlg") = g"^^^+'^'/' n (1-?''^'''^)"', (2.7) 

{i.j)e\ 
where (i,j) stands for the (i, j)-th box in the Young diagram, and n{X) is given 

by 



KA) = E(* - 1)^^ 



2.4 Melting crystal model with external potentials 

Wc now deform the foregoing melting crystal model by introducing the external 
potentials 



<|,,(A,p) = E 9'='^+^'-'+'^ - E '^"^"'^'^ 



1=1 i=l 

with coupling constants ifc, fc = 1, 2, 3, . . ., on the main diagonal slice A = 7r(0). 
The right hand side of the definition of $/c(A,p) is understood to be a finite sum 
(hence a rational function of q) by cancellation of terms between the two sums: 

cf>,(A,p) ^ ^(5'=b+>.-+i) _ gMP-^+D) + q^t^^. 
i=l ^ 

The partition function of the deformed model reads 

Zp(i)=E«'"'^*^*'"^°^'''^' (2.8) 

where 

$(i,A,p) = E^fe*fe(^'P)- 

fc=i 



We can repeat the previous calculations in this setting to rewrite the new par- 
tition function Zp(t) as 

Z,(t)=^.,(q^)2g*(*-^'^). (2.9) 

A 

Modifying this partition function slightly, we obtain the instanton sum of 
5D TV = 1 SUSY [/(I) Yang-Mills theory [6-. 

Zp{t) = J2 gl"lg"(°)e*(*''^(°)'P) = Y, SA(g'')'gl^le*(*'^'P). (2.10) 

q and Q are related to physical parameters R, A, H of 5D Yang-Mills theory as 

q = e-«^ Q = (i?A)2. 

The external potentials represent the contribution of Wilson loops along the fifth 
dimension |2Dj. In this sense, Zp{t)/Zp{0) is a generating function of correlation 
functions of those Wilson loop operators. 

Our goal is to show that the partition function Zp{t) is, up to a simple factor, 
the tau function of the (ID) Toda hierarchy. To this end, we now consider a 
fermionic representation of this partition function. 



3 Fermionic formula of partition function 

3.1 Complex fermion system 

Let i^iz) and ip*{z) denote complex 2D fermion fields 



oo oo 

r?i— — oo m — — oo 



4>{Z)= J2 ^rnZ-"'-', r{z)= J2 ^™ 



The Fourier modes ipm and ■(/'to of ip{z) and ?A*(z) satisfy the anti-commutation 
relations 

{V'm, CI = Sm+7i,0, {V'm, Ipn} = {V'm: O = 0- 

The Fock space F splits into charge p subspaces Fp-. 



F- e Fp. 



p— — oo 



The charge p subspace Fp has a unique normalized ground state (charge p 
vacuum) \p) and an orthonormal basis |A;p) labeled by partitions A. \p) is 
characterized by the vacuum condition 

V'mb) = for m > —p, 4'*n\p) — ^ ior m > p + 1. 



If the partition is of the form A = (Ai, . . . , A„, 0, 0, . . .), the associated element 
\X;p) of the basis is obtained from \p) by the action of fcrmion operators as 

|A;p) = V'-(p+Ai-l)-l • • • V'-(p+A„-n)-lV'(p-n) + l ' " ' i'{p-l) + l\p) ■ 

They are orthonormal in the sense that their inner products have the normahzed 
values 

(A;p|Ai;g) ^ SpgSx^. 

3.2 U{1) current and fermionic representation of tau func- 
tion 

The [/(I) current J{z) of the complex fermion system is defined as 

CXI 

k— — oo 

where : : denotes the normal ordering with respect to the vacuum 0) : 

:ipm1pn- = "iPmlPn " (0|V'mV'^ |0) . 

The Fourier modes 

'^m — / ^ •V^m — nV^n' 
n— — oo 

of J{z) satisfy the commutation relations 

[Jm, Jri] — ITl^m+nfi (3-1) 

of the Aoo Heisenberg algebra, and play the role of "Hamiltonians" in the usual 
fermionic formula of the KP and 2D Toda hierarchies [151 US] • For the case of 
tau functions T{t,t), t = (^1,^2, ■ ■ ■), t — (^1,^2, ■ • ■): of the 2D Toda hierarchy, 
the fermionic formula reads 

00 00 

Tp{t,t) = (p|exp(^ tmJni)gexp{- ^ tmJ-m)\p) (3.2) 

m— 1 771—1 

where g is an element of the infinite dimensional Clifford group GL{oo). 

3.3 Fermionic representation of Zp(t) 

The partition function Zp{t) of the deformed melting crystal has a fermionic 
representation of the form 

Zp(t)-(p|G+e^(*)G_b). (3.3) 



Let us explain the constituents of this formula along with an outline of the 
derivation of this formula. 

H{t) is the linear combination 



H{t)^J2^kHk 



of the "Hamiltonians" 



Hk= Y. '?'"-V'-nC:- 



n— — oo 



The aforementioned basis elements \X;p) of the Fermion Fock space turn out to 
be eigenvectors of these Hamiltonians. The eigenvalues are nothing but the the 
potential functions ^k{X,p)- 

Hk\X;p)^M\p)\^:p)- (3.4) 

G± are GL{oo) elements of the special form 

qk/2 



G±-e-p(Efc7T^7ry^±'= 



k=l 



Since the numerical factors q^/"^ /{\ — q^) in this definition can be expanded as 

_9 _ V^ -fc(m+l/2) _ V^ fc(m+l/2) 

I - qk Z^ ^ Z^^ ' 

m— — oo m— 

one can factorize these operators as 

— 1 oo 

G+= n r+M, G_-nr-M> (3-5) 

m— — oo m— 

where 

r±M=exp(^igT'=(™+V2)j^^y 
fc=i 

These T± (m) 's are a specialization of the so called vertex operators 

oo ^ 

V±{z) =exp(^^yJ±i 



for bosonization of the complex fermions. 

Following the idea of Okounkov and Reshetikhin [2] , we now consider m as 
a fictitious "time" variable. A plane partition then may be thought of as the 



10 



"path" (or "world volume") of discrete time evolutions of a partition A that 
starts from the empty partition = (0, 0, . . .) at infinite past and ends again in 
at infinite future (see Figure^]). The vertex operators T±{m) play the role of 
transfer matrices between neighboring diagonal slices. 

The vertex operators T±{m) act on the aforementioned orthonormal bases 
\X;p) and (A;p| as 

(A;p|r+(m) = ^(Ai;p|g-("+i/2)(|/.|-|A|) (3 g^ 

for m — —1, —2, . . . and 

r_(m)|A;p) - Yl -Z^^+'/'^^'^'^'^'V;?) (3.7) 

tj,yx 

for m — 0, 1, . . . [11 [3]. The right hand side of these formulas give a linear 
combination of all possible time evolutions of the rn-th slice A — 7r(m) at the 
next time. The weight gT(m+i/2)(|/j|-|A|) qJ each state on the right hand side 
are exactly the factors assigned to the boxes of 7r(m)/7r(TO =F 1) in the definition 
of the weights q^ , q^ that appear in the combinatorial formula (|2.4p . 

Since G± are products of these slicc-to-slice "transfer matrices", {p\G+ and 
G-\p) become linear combinations of the states {X;p\ and |A;p) that evolve from 
the ground states {p\ and \p) at m = ^cxd. By what we have seen above, the 
weights of {X',p\ and \X;p) in these linear combinations are given by the partial 
sums of q^ and q-^ over all semi-standard tableaux T and T' of shape A, namely, 
the special value s\{qP) of the Schur function. Thus {p\G+ and G-\p) can be 
expressed as 

{p\G+ = j2 E 'z^(A;pi-E^^(«')<^;pi' (3.8) 

A T:shapcA A 

G-\P)=T. E 'z'''|A;p)=5],s,(gP)|A;p). (3.9) 

A T':shapc A A 

The expectation value of e^^*-' with respect to these states yields the fermionic 
representation p.3p of the partition function Zp(t). 

The fermionic representation (j3.3|) is apparently different from the fermionic 
formula p.2[) of tau functions of the 2D Toda hierarchy. To show that Zp{t) is 
indeed a tau function, we have to rewrite p.3|) to the form of p.2p . This is the 
place where the quantum torus Lie algebra joins the game. 
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4 Quantum torus Lie algebra 

4.1 Fermionic realization of quantum torus Lie algebra 

Let Vm {k — 0,1, ... , m G Z) denote the following ferniion bilinear forms: 

n— — oo 

Note that 

•Jm — Vm , J^k — yQ 

Actually, Vm coincides with Okounkov and Pandharipande's operator £m{z) 
[5TJ [21] specialized io z — q^. As they found for f„i(z), our Vm s satisfy the 
commutation relations 

kA-l 

This is a (central extension of) q-deformation of the Poisson algebra of functions 
on a 2-torus. We refer to this Lie algebra as "quantum torus Lie algebra" . More 
precisely, a full quantum torus Lie algebra should contain elements for fc < as 
well; for several reasons, we shall not include those elements. 

4.2 Shift symmetry among basis of quantum torus Lie al- 
gebras 

The following relations, which we call "shift symmetry" , play a central role in 
identifying Zp{t) as a tau function: 



G_G+(l^« -<5„,0Y^)(G-G+)-i = i-lf{v!^l-5m+k.,^) (4.2) 



These relations are derived as follows. 

Let us recall that the fcrmion fields ^{z),^*{z) transform under adjoint 
action by J±fc's as 

expf^CfcJ±fcj-0(2)expf-^CfcJ±fej = eyL-pV^CkZ^^jipiz), (4.3) 

k=l k=l k=l 

oo oo oo 

exp(^CfcJ±fc)v*(2)exp(-^CfcJ±fc) = exp(^- ^ CfeZ±'=) V*(2)- (4.4) 



fc=i fe=i fe=i 
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By letting c^ = (^1"^ l{\ — q'^), the exponential operators in these formulas turn 
into G± , so that we have the operator identities 

G+^p{z)G+-' = (q^'^z;q)-Ji,{z), (4.5) 

G+r{z)G+-' = {q^'^z-q)^r{z). (4.6) 

G_V(z)G_-l = {q^'^z-^-q)-^i,{z), (4.7) 

G_V*(z)G_-l = {q'/^z-^-q)^r{z), (4.8) 

where (z;g)oo denotes the standard g- factorial symbol 

oo 
n=0 

We use these operator identities to derive transformation of the fermion 
bilinear forms :ip{q'''^z)ip*{q~'^''^z): under conjugation by G±. Sinee 

k/2 

■.i>{q^'^z)i>*{q-^'^z): = -i>* [q-^'^ z)i>{q^/^ z) + '^ 



(l-g'=)z' 

let us first consider ip* [q^^^"^ z)'ip{q^/'^ z) . Under conjugation by G+, it trans- 
forms as 



G+r{q-'''^z)^{q'''^z)G+-^ 
(qi/2 . q-^/'^z\q) 



(qi/2 . qk/2^.^q-j^ 



i;*{q-^/^z)i>{q'''^z) 



k 



= [] (1 - zg('=+i)/2-™)^*(g-'=/2z)^(gfc/2z). 

m— 1 

This implies that 

G^{-Mq'/'z)r{q-'/'z): ^J^)g+-^ 

= n(l-g('=+i)/^-"z)(:^(g'=/^z)r(g-'=/'^):-7Y^^). (4.9) 

In much the same way, we can derive a similar transformation under conjugation 
by G_. In this case, it is more convenient to rewrite the result as follows: 

G^-^{-Mq'/^zW{q-'/'zy.-^^^)G^ 

= n (1 - q-^'+'^/'+^z-'){-Mq''^'z)r{q-'/'zy. -^^)(4.10) 
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We note here that the prefactors on the right hand side of the last two 
equations are related as 

k k 

n (1 - 5(fe+i)/2-™^) = (_;^)fc [] (1 _ 9-(fc+i)/2+™^-i). 

m— 1 m—1 

Accounting for this simple, but significant relation, we can derive the identity 

= (-zf(:^(g^/^z)r(g~^/^z):- ^^^'^J,^J . (4.11) 

The shift symmetry (J4.2I) follows immediately from this identity. 
When m — and m ~ —k, (14.21) takes the particular form 

G-G+lv^"^ _^)(G_G+)-i = {-lM''\ (4.12) 

{G^G+r'lv,^"^ ^^)g^G+ = (-l)^Vit.^ (4.13) 

It is these identities that we shall use to convert the fermionic representation of 
Zp{t) to the standard fermionic formula of tau functions. 

5 Integrable structure of melting crystal model 

5.1 Partition function as tau function of 2D Toda hierar- 
chy 

Let us split the operator G+e^^*-'G_ in p.3p into three pieces as 

G+e^'^'^G- = G+e^(*)/2e^(*)/2G_ 

= G+e^(*)/2G+-i • G+G- ■ G_-ie^(*)/2G_ 

and use the special cases (|4.12p and ()4.13p of the shift symmetry to rewrite those 
pieces. 

To this end, it is convenient to rewrite (|4.12p and (|4.13p as 

k 

G+(fffc-^^)G+-i = (-l)''^G_-iyf)G_, 

k 

G--^[Hk-j^)G- = {-If G +¥["",! G+'\ 

Though the operators V^^, on the right hand side are unfamiliar in the theory 
of integrable hierarchies, we can convert them to the familiar "Hamiltonians" 
J-l-j, = V^jI of the Toda hierarchies as 

14 



where W is a special element of Woo algebra: 



n— — oo 



We thus eventually obtain the relations 



^A; 



G+[h,-^)g+-' = {-l)'G--\-'^/'j,q'^/'G-, (5.2) 

G_-i(i7,-^)G_ = (-l)'=G+g'^/V_fe<z-^/2G+-i (5.3) 

between HkS and J±kS. 

By these relations, G+e^*^*^/^G+^^ can be calculated as 

G+e^(*)/2G+-i 

OO ^ OO / ^ -.f. 

= -p(E27rr7^)G--.-'^/^exp(E(4i^J.),-/^G_. 



V^2(l-,'^). 



fe=i 



A similar expression can be derived for G_ ^e^*^*^G_ as well. We can thus 
rewrite G+e^'*^G_ as 



G,e-(*)G. =exp(^^)G_-,--/^exp(5:^^J.) x 

oo / ^ \ j^ , 

X5exp(E^^'/-.)'Z-^'/^G+-i 



(-l)'^ife 
fc=i 
where 

g ^ g^/2(G_G+)\'^/2. (5.4) 

The partition function Zp{t) is given by the expectation value of this operator 
with respect to {p\ and \p). Since the action by the leftmost and rightmost pieces 
of g yields only a scalar multiplier to (p|, \p) as 

(p|G_-ig-'^/2 = 5-p(p+i)(2p+i)/i2(p|^ (5.5) 

<Z-^/2G+-1|p) = g-p(p+i)(2p+i)/i2|p)^ (5.6) 

Zp(^) can be expressed as 

oo j^ 

Z,{t) ^ exp(^i^),-^(^+i)(^^+i)/«x 



fe=l ^ 



(p|exp(E^^^^)5exp(E^^^-^)|p)- (5-7) 



fe=i fc=i 



The expectation value (p| • • • \p) takes exactly the form of p.2p . Thus, up to the 
simple prefactor, Zp(i) is essentially a tau function of the 2D Toda hierarchy. 
Thus we find that an integrable structure behind the melting crystal model is 
the 2D Toda hierarchy. This is, however, not the end of the story. 
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5.2 ID Toda hierarchy as true integrable structure 

The foregoing calculation is based on the splitting 

G+e«(*)G_ = G+e^(*)/2G+-i • G+G^ ■ G^-ie^W/^G^. 

Actually, we could have started from a different splitting of G+e^'^*-'G_, e.g., 

G+e^(*)G_ = G+e^(*)G+"^ • G+G- = G+G- ■ G-^^e^^^^G- 

This leads to another set of expressions of Zp{t) in which only the {p\- ■ ■ \p) part 
is different. We thus have the following three different expressions for this part: 



00 

= (p|exp(E(-l)^-tfcJfc)5|p) 

fe=i 

00 

= {p\ge^^(Y,{-l)HkJ-u)\p). (5.8) 



fc=i 

These identities of the expectation values can be directly derived from the 
operator identities 

Jfe5 = 5J_fc, fc = 1,2,3,... (5.9) 

satisfied by g. (These operator identities themselves are a consequences of the 
shift symmetry of Vm s.) Generally speaking, this kind of operator identities 
imply symmetry constrains on the tau functions |231 124] : in the present case, 
the constraints read 

d , -, d 



—rp{t,i) + —Tp{t,t)^0, fc = 1,2,3,.... (5.10) 



In other words, the tau function is a function of i — i, 

Tp{t, t) ^Tp{t- t, 0) ^ Tp{Q, i - t), 

and reduces to a tau function Tp{t) of the ID Toda hierarchy that has a single 
series of time variables t — (ii, i2, ■ • •) rather than the two series of the 2D Toda 
hierarchy. Thus the ID Toda hierarchy eventually turns out to be an underlying 
integrable structure of the deformed melting crystal model. 

The same conclusion can be derived for the instanton sum (|2.10l) of 5D SUSY 
U{\) Yang-Mills theory. It has a fermionic representation of the form 

Zp(t) = (p|G+g^«e^«G_|p) (5.11) 

where Lq is a special element of the Virasoro algebra: 

00 

-^0 = X! ^■^-n'4'*n- 
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One can repeat almost the same calculations as the previous case to convert 
Zp{t) to the form of (|5.7p . The counterpart of g is given by 

g = q^/^G-G+Q^°G-G+q^^^, (5.12) 



which, too, satisfy the reduction conditions (|5.9p to the ID Toda hierarchy. 
Thus a relevant integrable structure is again the ID Toda hierarchy. 

6 Concluding remarks 

6.1 Problems on 4D instanton sum 

In deriving the instanton sum (|2.10p . 5D space-time is partially compactified in 
the fifth dimension as R'' x S^. The parameter R is the radius of S^. Therefore, 
letting i? — i> amounts to 4D limit. 

Unfortunately, it is not straightforward to achieve such a 4D limit in the 
present setup. Firstly, the 5D instanton sum with external potentials does not 
have a reasonable limit as i? — > 0. Any naive prescription letting R —^ yields 
a result in which t dependence disappears or becomes trivial |12| . Secondly, the 
shift symmetry of the quantum torus Lie algebra ceases to exist in the limit as 
q — e^^^ — )■ 1. Speaking more precisely, the quantum torus Lie algebra itself 
turns into a Woo algebra in this limit, but no analogue of shift symmetry (14.21) 
is known for the latter case. For these reasons, the 4D case has to be studied 
independently. 

The 4D instanton sum [71 18], too, is a sum over partitions. Moreover, this 
statistical sum has a fermionic representation [Ml [9] . Marshakov and Nekrasov 
[TTI [T8] further introduced external potentials therein. Actually, the 4D instan- 
ton sum for U{\) gauge theory is almost identical to the generating function of 
Gromov-Witten invariants of CP^ [HH^^]. This can be most clearly seen in the 
fermionic representation of these generating functions, which reads 

Zf{t) = (p|e^^/''exp(f;i,^)e'^-/''|p), (6.1) 

fc=i 

where Vk^ are fermion bilinear forms introduced by Okounkov and Pandhari- 
pande for a fermionic representation of (absolute) Gromov-Witten invariants 
of CP^ [21]. As regards these Gromov-Witten invariants (in other words, cor- 
relation functions of the topological a model) , it has been known for years 
[25l|26l [27l [28l|29] that a relevant integrable structure is the ID Toda hierarchy. 
Thus the ID Toda hierarchy is expected to be the integrable structure of 
the 4D instanton sum as well. This has been confirmed by Marshakov and 
Nekrasov in detail [T71 |T5]. What is still missing, however, is a formula like 
(|5.7p that directly connects Zf^{t) with the standard fermionic formula (13. 2p 
of the tau function. Finding a 4D analogue of (|5.7p is thus an intriguing open 
problem. This issue is also closely related to the fate of shift symmetry (|4.2p in 
the q — >■ 1 limit. 
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6.2 Relation to topological strings 

Our results are directly or indirectly connected with some aspects of topological 
strings as well. 

1. According to the theory of topological vertex [S], the partition function 
Zp{t) of the deformed melting crystal model has another interpretation as the 
A-model topological string amplitude for the toric Calabi-Yau threefold O © 
0{—2) — >■ CP"*". In this interpretation, q and Q are parametrized by the string 
coupling constant gst and the Kahler volume a of CP^ as 

Specializing the value of t leads to several interesting observations [T2] . 

2. A generating function of the two-legged topological vertex W\fj, ^ c\^, is 
known to give a tau function of the 2D Toda hierarchy [3D]. In the fermionic 
representation p.2p . this amounts to the case where 



9 



= q'^'^G+G.q'^'\ (6.2) 



(Actually, for complete agreement with the usual convention, we have to replace 
W with 

1\2 



^= E ("-2) ^'^-"^-' 



but this is not a serious problem. The difference can be absorbed by rescaling 
ifc's.) Let us stress that this GL(oo) element does not satisfy the reduction 
condition (|5.9I) to the ID Toda hierarchy. 

3. A generating function of double Hurwitz numbers for coverings of CP^ 
gives yet another type of tau function of the 2D Toda hierarchy [31] . Actually, 
the GL{oo) element for the fermionic representation is given by 

g = q^/\ (6.3) 

(More precisely, as in the previous case, W has to be replaced by K, but the 
difference is again irrelevant.) In this case, the reduction condition ()5.9|) to the 
ID Toda hierarchy is not satisfied, but the operator identities (15. ip imply that 
another set of reduction conditions are hidden behind (see below). 

6.3 Constraints and quantum torus Lie algebra 

As a consequence of the shift symmetry of Vm 's, the GL{oo) elements g of the 
aforementioned models of topological strings turn out to satisfy some algebraic 
relations other than (|5.9p . According to general results on constraints of the 2D 
Toda hierarchy [231 HI] , such relations imply the existence of constraints on the 
tau functions and the Lax and Orlov-Schulman operators. Those constraints 
inherit the structure of the quantum torus Lie algebra. Let us illustrate this 
observation for the case of double Hurwitz numbers over CP^. 
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The GL{oo) element g — q^l"^ for this case satisfies the operator identities 

^fc5 = 5^i'\ 9J-k = y^^^9 (6.4) 

as a consequence of (|5.1I) . These identities can be converted to the constraints 

L = q^'\^'L, Z-1 = g-i/V'i-1 (6.5) 

on the Lax and Orlov-Schulman operators L, M, L, M of the 2D Toda hierarchy. 
Emergence of the exponential operators q^ and q^ is a manifestation of the 
quantum torus Lie algebra. To see this, let us recall that the Lax and Orlov- 
Schulman operators satisfy the (twisted) canonical commutation relations 

[L, M] = L, [L, M] = L. (6.6) 

This implies that the monomials q-k^n/2j^,nqkM ^^^^ ^-km/2im^kM ^f ^^ ^m ^ i^ ^m 
give two copies of realizations of the quantum torus Lie algebra. 

The constraints (j6.5p are remarkably similar to the "string equations" 

L = ML, L^^ ^ ML^^ (6.7) 



of c = 1 strings at self-dual radius [32l [33l [23l [24] . A relevant algebraic structure 
of these string equations is the Woo algebra. Thus (|6.5p may be thought of as 
g-deformations of these M^-algebraic constraints. 
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